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Context. Long-period variables (LPVs) and Miras exhibit a pronounced variability in their luminosity with a more or less well-defined 
period, and they suffer large mass loss in the form of stellar winds. Owing to this extensive mass loss, they are surrounded by extended 
circumstellar dust shells (CDSs). The dynamics of these envelopes is the result of a complex interplay via an external excitation by 
the pulsating central star, dust formation, and radiative transfer. 

Aims. Our study is aimed at an understanding of the dynamics of CDSs around carbon-rich, standard luminous LPVs and Miras. 
These shells often show multiperiodicity with secondary periods as high as a few 10 4 d superimposed on a main period that is in the 
range of approximately 10 2 - 10 3 d. Such secondary periods may be caused at least in part by the presence of dust. 
Methods. We consider an excitation of the CDSs either by a harmonic force, provided by the oscillation of the central star, or by a 
stochastic force with a continuous power spectrum. The resulting numerically computed dynamical behaviour of the shell is analysed 
with the help of Fourier analysis and stroboscopic maps. 

Results. CDSs may be described as multioscillatory systems that are driven by the pulsating stars. A set of normal modes can be 
identified. The obtained periods of these modes are some 10 3 d, which is a characteristic timescale for dust nucleation, growth, and 
elemental enrichment in the dust formation zone. Depending on the oscillation period and strength of the central star, the envelope 
reacts periodically, multi- periodically, or irregularly. 
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1. Introduction 

Long-period variables (LPVs) and Miras are highly evolved 
stars located in the upper region of the asymptotic giant branch 
(AGB). They are characterised by low surface temperatures T < 
3 x 10 3 K and high luminosities of the order of 10 4 L . They typi- 
cally exhibit a pronounced variability of their intensity and suffer 
large mass loss in the form of a stellar wind with mass loss rates 
of up to 10~ 4 M G yr~' . Owing to their extensive mass loss, LPVs 
and Miras are surrounded by extended circumstellar envelopes 
that provide the necessary conditions for the formation of dust. 

|Dreyer et al. ((2009') (henceforth called paper I), a previous 
paper of this series, studied the dynamics of circumstellar dust 
shells (CDSs) around high luminous LPVs and Miras. It turned 
out that for these kind of stars, the dynamics of the envelopes is 
clearly dominated by the interplay between dust formation and 
radiative transfer, i.e. the exterior K-effect (see Flei scher et al. 
|1995||Hbfner et al.| 1995 i. A regular oscillation of the shell can 
be solely explained by the formation of dust and the radiative 
pressure on it. 

In this work we concentrate on the circumstellar envelopes 
of less luminous LPVs and Miras. For luminosities of less than 
10 4 Lq, additional input of energy and momentum is required to 
sustain the observed oscillations. This energy is provided by os- 
cillations of the star's deeper layers. Caused by the pulsational 
instability due to variations in the convective flux (cf. |Fox & 
|Wood |(T982[ l, hydrodynamic waves are generated in the dense 
stellar atmosphere and steepen into shock waves when travel- 



ling into the thinner circumstellar medium. The presence of such 
shock waves causes a levitation of the atmosphere producing 
favourable conditions for dust nucleation and growth. Once the 
first dust particles have formed, the momentum coupling be- 
tween matter and radiation field is strongly enhanced, which 
again results in a massive stellar wind. 

Although the shell's oscillation has to be triggered by a stel- 
lar pulsation, the presence of dust influences the dynamics of the 
entire system. Even for a strictly harmonic excitation, the shell 
does not necessarily mimic the stellar pulsation, but instead ex- 
hibits its own normal modes (cf. Bar dies & Tuchman ||1994) . 
Superimposed upon the excitation period that is in the range of 



approximately (10 2 - 10 3 ) d (jwhitelock et al. 



1991 



, the envelope 

typically shows oscillations with periods of up to a few 10 4 d. 
This could explain, at least in part, the observed cycle-to-cycle 
variation in the periods of some Miras fEddington & Plakidis~| 
T929}. As pointed out by|Percy & Colivas |(|1999[), these fluctua- 
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tions could come from multiperiodicity, with either a secondary 
period that is comparable to the primary one or a secondary pe- 
riod an order of magnitude longer than the primary one. In the 
latter case, the secondary oscillation may be explained by an os- 
cillation of the circumstellar shell from the exterior K-effect. 

In the context of dynamical systems, the CDS can be de- 
scribed as a multioscillatory system that is driven by the pulsat- 
ing stellar atmosphere. The resulting dynamical behaviour de- 
pends on both the frequency and strength of the excitation, as 
well as, on the intrinsic dynamics of the shell. For a close ex- 
amination of the possible normal modes of a CDS, which are 
controlled by the intrinsic timescales of the various physical and 
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chemical processes, we consider response spectra of a circum- 
stellar shell that is excited not only by a harmonic force but 
also by a stochastic force with a continuous power spectrum. 
The inherently non-linear nature of the underlying equations in- 
troduces an element of randomness into the shell's dynamics. 
This also becomes clear when the CDS is excited by a strictly 
harmonic force. For a sensitive range of frequencies and excita- 
tion strengths, the shell does not behave strictly periodically but 
shows signs of deterministic chaos. 

This paper is organised as follows. Section [2] outlines our 
dynamical modelling method of a CDS, giving a special focus to 
the formulation of a purely harmonic and stochastic excitation of 
the CDS provided by a central star. Section[3]presents our direct 
methods as used to investigate the dynamical CDS behaviour. 
This section discusses the problems of creating and analysing 
discrete power spectra and deterministic maps. In Sect. [4] we 
apply the method to a reference model. The eigendynamics 
are first determined and then the complex interaction between 
the inner excitation by a pulsating star and the multioscillatory 
circumstellar shell is discussed, in this case depending on 
excitation frequency and strength. The paper finishes with some 
concluding remarks in Sect. [5] 



2.2. Stochastic excitation 

To examine the eigenfrequencies and the band-pass character- 
istics of circumstellar shells around standard and low-luminous 
AGB stars, we consider the excitation of the shell by external 
white noise as generated by, e.g., a convective atmosphere of the 
giants. In analogy to the piston approximation, the movement of 
the inner boundary /?; n (f) is described by the following stochastic 
equations: 
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Rin(t)=R + r(f')df' 



M in (?) = r(o 



(3) 



(4) 



The term F(f) describes the influence of the randomness or noise 
and is taken to be a Gaussian with a mean of zero: 



<r> = o. 



(5) 



Since we adopt white noise, T(f) is uncorrected in time, which 
corresponds to the (auto)correlation function K(t) being Dirac's 
5-function 



2. Modelling of CDSs 

The theoretical model of a CDS is conceived as the complete 
unique set of solutions of an appropriate system of fundamen- 
tal dynamical equations, comprising the conservation laws, 
transport equations, and constitutive equations that describe 
both the local and global physical situation of the CDS and the 



various inherent processes, such as radiative transfer (Lucy 
TWTI [T9761 |Unno & Kondo ||1976|l or dust formation ( |Gail & 



Sedlmayr 



1988; 



Gauger et al. Q990). These equations were 



listed and solved according to the procedure given in |Fleischer] 



et al. 



T992J. The CDSs are characterised by the fundamental 
parameters of their generating stars, which are mass M„, 
luminosity L*, effective temperature T*, and carbon-to-oxygen 
abundance ratio C/O. As these parameters change only slightly 
with time because of stellar evolution, they are considered as 
given external parameters for modelling a circumstellar shell. 



K(t) = (T(t)Y(t + r)> = ct 2 6(t) , 



(6) 



where cr 2 scales the intensity of the noise. According to the 
Wiener-Khinchin theorem, the spectral density \H(f)\ 2 of T(f) 
is given by the Fourier transform of its autocorrelation function, 
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and is therefore independent of the frequency /. The integral 
of the Gaussian white noise, i.e. the solution of Eq. Q, is 
known as the Wiener process. It is also Gaussian distributed, 
so the mean vanishes, and its variance increases linearly with 
time t. For a concrete numerical modelling of the CDS, only 
realisations of the stochastic process Eq. Q where R m {t) stays 
within predefined finite boundaries (to ensure that the inner 
boundary stays within the stellar radius R t ) are considered. The 
strength of the noise cr is also an input parameter. 



2.1. Monoperiodic excitation 

The pulsation of an AGB star, more specifically the underly- 
ing inner /e-mechanism, is treated by the piston approximation 
adopted from the work of|Wood | ( |l979| l, |Bowen | ( |1988[ ), and 
|Bowen & Willson~| ( |1991| l. It is assumed that the radius of the 
inner boundary Ri n (t) located a few scale heights below the stel- 
lar radius oscillates sinusoidally about an equilibrium position 
/?o with 



P l2n 
Rinit) =R +Au— sin \—t 



(1) 



This implies that, at the inner boundary the hydrodynamic ve- 
locity M; n , varies in a manner determined by 



>2n 

u m (t) = Am cos I —t 



(2) 



The velocity amplitude Am, which simulates the strength of the 
interior pulsation of the star, as well as the pulsation period P, 
are input parameters. 



3. Diagnostic approach 

The dynamics of (non-)linear systems may be described by 
means of different characteristics. One of the classical ways to 
characterise the dynamics is to compute the power spectrum 
based on the Fourier transform. This provides the possibility of 
obtaining a deeper physical understanding of the real system dy- 
namics by detecting various frequencies contained in a given sig- 
nal. Another way is an examination by means of a deterministic 
approach, e.g., to create stroboscopic or specific Poincare maps 
to find domains of stability in the phase space. 

3. 1 . Fourier transform 

Generally, the Fourier transform is simply a method of express- 
ing a function in terms of the sum of its projections onto a set of 
basis functions. The standard basis functions used for Fourier 
transform are {sin(27r/f), cos(2nft),f e R) or, equivalently, 
{e 2m f',f G R}. It is the oscillatory frequency / that varies over 
the set of all real numbers to give an infinite collection of basis 
functions. We project a given time function h(t) onto our basis 
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functions to get the Fourier amplitudes H( f) for each frequency 
/: 



(8) 



In general, H(f) will be complex. The norm of the amplitude, 
\H(f)\ is called the Fourier spectrum, and the square \H(f)\ 2 is 
called the power spectrum of h(t). 

When using a fast Fourier-transform (FFT)-algorithm to 
compute the spectrum of a signal, we use a finite length of sig- 
nal, consisting of N = 2 y ,y e N discrete samples. However, the 
discretisation of a signal is not without problems. If the sampling 
rate / s = 1 /At is not high enough to sample the signal correctly, 
then a phenomenon called aliasing occurs; i.e., components of 
the signal at high frequencies are mistaken for components at 
lower frequencies. 



3.2. Mapping 

For dynamic systems having several degrees of freedom, it is not 
very practical to discuss the orbit in a multidimensional phase 
space. A more appropriate way is to study the intersections of 
orbit with a plane in phase space, which reduces the phase space 
dimension by one. 

If the differential equations describing a system contain 
terms periodic in time (where P is the period), it is convenient to 
create stroboscopic maps showing the time trajectory for several 
time intervals Af; i.e., the flow is sampled whenever t - nAt for 
n e Nq. This kind of map can be constructed for any temporal 
interval. Sampling by a fixed time interval Af = P, i.e. by suc- 
cessive times when the orbit crosses the phase angle <p plane, i.e. 
t = kin with k being an integer, results in a Poincare map. 



set of fundamental stellar parameters, but for different realisa- 
tions of the stochastic process as described in Sect. 2.2 The 



stochastic movement of the inner boundary R m {t) and the cor- 
responding velocity u\ n (t) for an arbitrary stochastic realisation 
with intensity cr = 0.1 kms -1 are displayed in Fig.[T] 




timet [d] 



Fig. 1. Stochastical movement of the inner boundary about an 
equilibrium position Rq = 0.87 R t (upper panel) and associated 
velocity u m (lower panel) for the reference model, stochastically 
excited with strength cr = 0.1 kms -1 . 



4. Application 

The condition for a purely dust-driven wind by the exterior k- 
mechanism is only properly fulfilled at the very end of the AGB 
evolution. If the stellar luminosity decreases, the dominance of 
the dust with respect to the CDS dynamics at first diminishes. 
Then, further reduction eventually leads to a vanishing out- 
flow, if no additional energy input is provided. AGBs with low 
and standard luminosities generate neither self-induced shocks 
and outflows, with self-maintained oscillatory patterns caused 
solely by the dust formation (exterior /(--mechanism), nor a CDS. 
Generating a CDS requires input of momentum and energy via 
an additional mechanism in order to levitate material into the 
outer atmosphere, i.e. into the dust formation window. An exam- 
ple of such a mechanism can be the pulsation of the outer layers 
of a red giant. The pulsation provides the star with mechanical 
momentum for levitating the material in their atmosphere. 

4. 1 . Eigendynamics of envelopes 

To examine the eigendynamics of a CDS around a standard or 
low luminous AGB star (e.g. with the reference stellar param- 
eters: L, = 9 x 10 3 L , M, = 1 M , T, = 2, 5 x 10 3 K, and 
C/O = 1.75), the additional mechanism of pulsation is adopted 
to be irregular (stochastically determined), such as is generated 
by convection in a red giant atmosphere. The reason for this as- 
sumption is to avoid the dominance of the stellar pulsation on the 
CDS dynamics. The shell has been excited with different noise 
intensities cr = (0.1,0.5, 1.0,2.0) kmr 1 and sequences from up 
to 24 models were calculated (cf. Table [TJ, each with the same 



Figure [2] shows snapshots of the radial envelope structure 
from the stochastically excited model calculation at time t = 
38 800 d for the intensity strengths cr = O.lkms -1 (left) and 
cr = l.Okms -1 (right) for various physical quantities describ- 
ing a CDS. The higher input of mechanical energy in the model 
with cr = l.Okms -1 can be seen immediately from the higher 
amplitude of the innermost shock compared to the model with 
cr = O.lkms -1 (uppermost panels). Since the density scale 
height is the same for both models in the inner region, the higher 
energy input yields a smoother transition to the wind region 
(second panels), i.e. a more levitated atmosphere, hence more 
favourable conditions for dust formation. The outflow is initi- 
ated by the pulsation but is accelerated by the dust. The stel- 
lar radiation is blocked by the opacity of the instantly formed 
shell, thereby heating the material behind it, which is located at 
smaller radii. This backwarming effect prevents new dust from 
forming and demonstrates that formation and growth of grains 
also strongly influences the thermodynamical shell structure. As 
depicted in the lowermost left panel for cr = 0.1 km s -1 , the ra- 
dial structure is rather smooth with an almost constant degree 
of condensation / conc j « 0.2; i.e., only 20% of the condensable 
carbon has actually condensed into solid particles. 

This results in an optically thin, slowly dust-driven wind 
with a time-averaged mass loss rate of (Mrp 9.5 x 10~ 8 M yr~' 
and terminal velocity of ~ 3.6 km s~ . For cr — 1.0 km s _1 , 
the wind becomes rather irregular, but shows the well-known 

1 (x) = l/T f t ' +T x(f') dt'. This nomenclature denotes averaging of 
quantity at a fixed radial position of r = 25 R, over a time interval T of 
at least one pulsation period P. 
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t = 38 800 d 



t = 38 800 d 





Fig. 2. Radial structure of stochastically excited CDSs with excitation strength cr = O.lkms -1 (left) and cr = l.Okms -1 (right). 
Upper panel: gas velocity u (solid red line), gas temperature T (dashed blue line). Second panel: mass density p (solid orange 
line), radiative acceleration on dust in units of the local gravitational deceleration a (dashed grey line). Third panel: stationary 
nucleation rate J* (solid light green line), number of dust grains per H-atom Hd/«H (dashed dark green line). Lower panel: degree of 
condensation of carbon / con( j (solid black line), mean particle radius (a m ) (dashed cyan line). 



onion-like shell structure (lowermost right panel). The calcu- 
lated mass loss rate is (M) = 4.5 x 10~ 6 M yr _1 , and the final 
outflow velocity shows a typical value of several tens of kilo- 
metres per second for an AGB star, namely (m m ) * 25.4 kms~'. 
For both excitation strengths the nucleation zone is localised 
around (2 - 3) R* (third panels). 

These stochastically excited CDSs were subjected to a 
Fourier transform. Since the gas velocity u(r, t) of a CDS outflow 
is usually taken to be a tracer of (multi)periodicity, we focused 
on a detailed analysis of the velocity structure. Regularly sam- 
pled with an interval At over a time domain of t max = 160 000 d, 
the time-dependent velocity u(r, t) was transformed at each fixed 
radial position r for 1 R„ < r < 5 R„, i.e. within the dust-forming 
zone. The spatial resolution was chosen to be Ar = 0.25 R*. To 
identify whether spurious frequencies are generated either by 
interaction between the grid with numbers of gridpoints j max , 
grid remapping rezoin<r\or by the choice of the sampling in- 
terval Af, we studied several numerical parameter combina- 
tions. These studies suggest the following optimal numerical 
parameter set for our reference CDS: N = 512, y max = 4098, 
rezoinc = 0.4 P, and At = 83 d, which yields a Nyquistfrequency 
of f c — 6 x 10~ 3 d~'. Averaging all stochastic realisations pro- 
vides an overall mean power spectrum of the local expansion 
gas velocity u, which is presented in Fig. [3] for the same ex- 
citation strengths cr = (0.1,1.0) kms~' as discussed in Fig. [2] 
See also Table [T] which lists the model parameters of excita- 
tion strength cr, the number of ensemble members No., and 
the resulting quantities: ensemble-averaged CDS-eigenmodes 
f K , time-averaged mass loss rate (M), outflow velocity («<»), 
and dust-to-gas mass ratio (pd/p g ) in the dust nucleation zone 
atr»2.5R H ,. 



For details of grid redistribution see Fleischer 



1994) 



To gain more insight, Fig. [4] illustrates the power spectra 
for some selected radial positions within the circumstellar en- 
velope. For small excitation strength cr = 0.1 kms~', the re- 
sponse spectra are dominated by three distinct peaks at f K - 
0.06 x 10 3 d 1 * (HOOOdr 1 , at f K = 0.51 x lO^d" 1 * 
(1960d)-\ and at f K = 1.58 x lO^d" 1 * ^SOd)" 1 . The am- 
plitude of the third frequency reaches its maximum near the 
inner edge of the dust nucleation zone at r/R t = 1.75 and di- 
minishes rapidly with increasing distance from the star. Beyond 
r/R t - 1.75, dust formation sets in and the first and sec- 
ond maxima at f K = 0.06 x lO^d" 1 * (17 000 d) -1 and at 
f K = 0.51 x 10~ 3 d l x (1 960 d)" 1 become visible. They eventu- 
ally dominate the spectrum for r/R t > 1.75. 

For the higher excitation strength cr = 1.0 kms~', the power 
spectrum looks much richer in detail than for cr = 0.1 kms 1 
with at least five distinct maxima. Even closer to the star where 
the medium is too hot to allow for effective dust formation, two 
strong peaks at /, = 1.13 x 10~ 3 d~' * (885 d)" 1 and f K = 1.50 x 
10~ 3 d~' * (670 d)~' can be seen. Here, the dynamics of the gas 
are influenced by the backwarming of the dust and the dilution 
wave originating in the dust nucleation zone (cf. Paper I). With 
increasing distance from the star, additional maxima appear at 
f K = 0.06 x lO^d- 1 * (17000d) -1 , f K = 2.13 x lO^d" 1 * 
(470 d)->, and/* = 3.2 x lO^d" 1 * OlOd)- 1 . 

Similar to the case of a self-excited CDS (Paper I), which 
can be characterised by an eigenmode P K , these frequency max- 
ima are characteristic of the dynamics of envelope around low 
and standard luminous LPVs and Miras. The result shows that 
a CDS does not simply mimic an external excitation but has its 
own eigenmodes, which are determined by the complex inter- 
play between dust formation, growth, radiative transfer, and hy- 
drodynamics. It has also become apparent that the response of 
the CDS strongly depends on the excitation strength. It is as- 
sumed that the maxima at f K = 1.58 x 10~ 3 d~' * (630 d) -1 
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Values in boldface indicate the dominant mode 




Fig. 3. Normalised mean power spectra of the radial gas velocity u caused by stochastically excitation of intensity <x = 0.1 kms 
(left) and cr = 1.0 km s" 1 (right). 



for cr = 0.1 kms" 1 and at f K = 1.50 x lO^d" 1 * (670 d)" 1 for 
cr = 1.0km s _1 belong to the same eigenmode, whose frequency 
is detuned towards lower values as the strength of the excitation 
increases. 



4.2. Dynamics of periodically excited envelopes 

This section examines the response of the CDS to an external 
excitation by a sinusoidally oscillating stellar atmosphere 
(cf. Sect. |2.1| i. Owing to the inherent non-linearities of a dust- 
forming system, the CDS does not generally oscillate with the 
applied frequency but usually exhibits a much more complex 
behaviour. Depending on the excitation frequency, it may 
exhibit an irregular or periodic behaviour. Since the dynamics in 
CDSs are strongly influenced by the presence of dust, we focus 
our investigation not only on the dust-formation zone, but also 
on the physical quantities responsible for providing favourable 
conditions for dust formation, such as gas velocity u, density p, 
and temperature T. 



4.2.1. Excitation frequency 

Figure [5] shows the most prominent shell modes normalised by 
the eigenfrequency fcos/fx m the dust nucleation zone at r « 
2.5 R, as a function of the excitation period P in the case of our 
reference model. Once again, the perturbation force was taken 
as the sine of period P according to Eqs. ([T} and Q, whereas the 
amplitude was held at Am = 1.0 kms -1 . 

The excitation frequency / p = l/P itself and some of its 
harmonics are recognisable for all periods P, whereas the deter- 
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Fig. 5. Most dominant frequencies of the CDS /cds normalised 
by the eigenmode f K = 0.51 ■ 10~ 3 d -1 (triangles), periodically 
excited with various stellar pulsation periods P but fixed ampli- 
tude Am = 1.0 kms -1 , at the dust nucleation zone r ~ 2.5 R*. 
The set of eigenfrequencies are indicated by circles. The eigen- 
period is labelled with P K . The excitation frequencies lie on the 
solid line and their harmonics on the dash-dotted lines. The grey 
line separates the irregular (left) from the pulsation-dominated 
domain (right). 



mined eigenmodes f K are only visible for small excitation peri- 
ods f<1250d. 
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Fig. 4. Selected mean power spectrum profiles of the gas velocity u at different radii r for the noise intensities cr -0.1 km s 1 (upper 
panels) and cr = l.Okms -1 (lower panels). 



- Irregular domain 

For small and intermediate excitation periods P < 1 250 d, 
the power spectra show a rather extended continuum with 
a number of distinct maxima. Some of these could be 
identified as eigenmodes of the CDS or as their harmonics. 
Also, the subharmonics f p /2 and f v /3 of the excitation 
frequency / p = l/P show up. As the pulsation period P 
increases, the excitation period and its harmonics become 
clearly recognisable. In contrast, the amplitude of the enve- 
lope's eigenmode mostly decreases as P increases. Within 
this domain a kind of resonance phenomenon between 
excitation and envelope period occurs. The envelope reacts 
multiperiodically for a sensitive range of excitation period P 
and strength Am. 

- Pulsation-dominated domain 

For longer periods P > 1 250 d, the CDS becomes finally 
enslaved by the excitation force. The most prominent feature 
of the spectrum is the pulsation mode supplemented by a set 
of its harmonics. 

In general, the final velocity {u^) has a fairly high value with 
a maximum at (Moo) = 32kms~ 1 for P - 600 d. With increasing 
period P, the velocity continuously slopes downward to (m m ) = 
19kms _1 at period P = 4500d. The corresponding mass loss 
rates show the same trend. The maximum mass loss rate is (M) 
1 .8 x 10~ 5 M yr 1 for P = 600 d, and the minimum (M) — 2.6 x 
10~ 6 M yr 1 for P = 4500d (cf. Table |}. 

Figure [6] demonstrates the basic behaviour of the enve- 
lope for some selected monoperiodic excitation periods P = 
376 d, 450 d, and 716 d chosen from the irregular shell dynam- 
ics regime. The left panels display the power spectra of the local 
expansion velocity u by scaling to the overall strongest mode. 
Finally, the right panels present the stroboscopic maps in the (u, 
p, T)-phase space, supplemented by a projection onto the (m, p)- 
plane. Both were sampled at an interval Af = 0.02 P. They also 



Table 2. List of resultant quantities from a CDS periodically dis- 
turbed. 
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9.4 x 10- 


-6 


25.6 


3.5 x 10" 


-3 


376 


1.33 


2.66 


3.99 


1.0 x 10" 


-5 


26.4 


4.1 x 10" 


-3 


450 


0.44 


0.65 


1.35 


1.1 x 10- 


-5 


25.8 


3.9 x 10" 


-3 


555 


0.39 


0.58 


1.19 


1.6 x 10" 


-5 


26.2 


4.1 x 10- 


-3 


600 


1.67 


3.33 


5.00 


1.8 x 10- 


-5 


31.8 


4.2 x 10- 


-3 


650 


0.77 


1.54 


3.08 


7.9 x 10 


-6 


28.9 


4.2 x 10- 


-3 


716 


1.40 


2.79 


4.19 


7.5 x 10- 


-6 


26.8 


4.2 x 10- 


-3 


900 


0.45 


0.54 


1.11 


9.7 x 10 


-6 


25.2 


3.2 x 10- 


-3 


1000 


0.35 


0.51 


1.00 


1.7 x 10- 


-5 


24.4 


2.7 x 10- 


-3 


1250 


0.54 


0.80 


1.60 


1.3 x 10- 


-5 


26.3 


4.2 x 10- 


-3 


1500 


0.67 


1.33 


2.00 


5.0 x 10" 


-6 


26.8 


3.7 x 10- 


-3 


1750 


0.57 


1.21 




4.3 x 10- 


-6 


25.7 


4.1 x 10- 


-3 


2000 


0.50 


1.00 


1.17 


3.8 x 10" 


-6 


23.5 


2.3 x 10- 


-3 


2300 


0.34 


0.87 


1.05 


2.1 x 10- 


-6 


23.1 


2.6 x 10- 


-3 


2500 


0.40 


0.80 


1.20 


1.8 x 10 


-6 


18.8 


1.4 x 10- 


-3 


2600 


0.38 


0.77 


0.89 


2.2 x 10- 


-6 


18.4 


1.9 x 10- 


-3 


2700 


0.37 


0.74 


1.09 


3.1 x 10- 


-6 


20.4 


2.9 x 10- 


-3 


3000 


0.33 


0.67 


1.01 


4.5 x 10- 


-6 


20.4 


2.7 x 10- 


-3 


3250 


0.31 


0.62 


0.92 


3.6 x 10" 


-6 


16.8 


2.3 x 10- 


-3 


3500 


0.27 


0.65 


0.85 


4.8 x 10- 


-6 


18.1 


1.9 x 10- 


-3 


4000 


0.25 




0.75 


2.6 x 10- 


-6 


18.8 


2.9 x 10- 


-3 



Values in boldface indicate the most dominant mode and underlined 
values the excitation frequency 



show the corresponding state of the system at a constant phase 
angle <p, i.e. At - P (Poincare map) in the dust nucleation zone 
atr « 2.5 R«. 



C. Dreyer et al.: Circumstellar dust shells around long-period variables X. 



7 





Fig. 6. Power spectra of the local expansion velocity u (left panels) and corresponding maps of the (u,p,T) phase space (green) 
and projection onto (M,p)-plane (yellow) (right panels) in the dust nucleation zone for different monoperiodic excitation peri- 
ods P = 716 d (upper), 450 d (middle), and 376 d (lower), respectively. The stroboscopic maps were obtained by sampling 
[u(n0.02P),p(n0.02P), T(n0.02P)] (green) and the Poincare maps by [u{nP), p(nP), T(nP)] (black) for 1 < n < [t max /P], n e N. The 
excitation frequencies are highlighted in the spectra. 



As can be seen in the uppermost and lowermost spectra, the 
CDS reacts with the same dominant mode / = 1 .34 x 10 3 d _1 = 
(716 d)~', even though it is excited with different periods P = 
716d and P = 376 d. This is similar to the resonance phe- 
nomenon (iP = v'Pcds) studied in paper I for high-luminous 



Miras and LPVs. However, for the case of standard luminous 
LPVs and Miras studied in this paper we always found 7=1, 
which means that the period of the CDS Pcds appears as an in- 
teger multiple of the excitation period P (iP = Pcds)- In the 
context of circumstellar envelopes, this phenomenon is often re- 
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ferred to as multiperiodicity (cf. |Winters et al.|| 1994} |Fleischer| 
et al. 1995 1 Hofner et al. 1996 1. As indicated by the number 
of the black clusters in the stroboscopic maps, the system re- 
acts monoperiodically (single black cluster, upper right panel, 
i.e. i = 1) for an excitation with P — 716 d, whereas the system's 
response is double-periodic (two black features in the lowermost 
plot on right, i.e. i = 2) for an excitation with P = 376 d. In other 
words, the system returns to the same state after a time interval 
Af = iP, The corresponding curves in the (u,p, T)-hyper plane 
are restricted and closed. Nevertheless, some parts of the trajec- 
tory are passed through extremely rapidly, so that these parts are 
only depicted by a few points on the corresponding stroboscopic 
map. 

For P = 450 d, as can be seen from the power spectra 
in the middle-left panel, the excitation period does not dom- 
inate the envelope. The amplitude of the piston frequency 
/ p = 2.22 x 1CT 3 d -1 is weak, and the frequency spectrum shows 
up with distinct maxima at / = 0.44 x 10~ 3 d _1 « (2300 d)" 1 , 
which is a characteristical timescale of the CDS and 
/ = 1.35 x 10~ 3 d~' * (740 d)" 1 , which seems to corre- 
spond to the somewhat-shifted feature P = 716 d providing the 
resonance case for this model. This is also consistent with the 
stroboscopic map given in the middle-right panel. In contrast 
to P = 716 d and 376 d, the system is not limited to a closed 
trajectory in the (u,p, T)-space but rather completely fills out a 
continuous part of the phase space. However, the dynamics of 
the envelope once more seem to synchronise with the excitation 
period to some extent. For a constant phase angle <p of the stellar 
pulsation, i.e. Af = P, the system stays inside well-defined 
stripes that cut through the entire phase-space area filled out by 
the system. 



4.2.2. Excitation strength 

This section studies the above-mentioned reference model with 
regard to the response of the CDS for a sinusoidally stellar 
pulsation with period P = 376 d, i.e. a typical oscillation pe- 
riod of LPVs and Miras, over a sequence of varying strengths, 
i.e. amplitudes Am. The amplitude starts with a value of Au = 
l.Okms -1 , so the disturbance is subsonic, and is then increased 
up to Au = 9.0 km s -1 i.e. into the supersonice regime. Table [5] 
lists the results of these calculations. 

The excitation frequency / p = 2.66 x 10 3 d" 1 = (376 d)" 1 
itself is only noticeable for a low excitation strength of Au = 
l.Okms -1 . Because weak in power, it vanishes with increasing 
excitation strengths. The inertial envelope is not able to follow 
the short excitation period; in fact, this mode at no time does 
it dominate the envelope dynamics. Actually, the first subhar- 
monic of the excitation frequency / p /2 = 1.33 x 10~ 3 d -1 w 
(752 d) -1 dominates the dynamics. It can be regarded as a de- 
tuned envelope-eigenmode f K = 1.50 x 10~ 3 d a (670 d) -1 . 
Subsequently, the eigenmode at f K = 0.50 x 10~ 3 d -1 w 
(2000 d)~' takes over the dominance of the shell dynamics. In 
general, the CDS-eigenmodes are detuned towards lower fre- 
quencies as already shown in the case of a CDS around a high- 
luminous LPV (cf. paper I). The stellar pulsation leads to a lev- 
itation of the atmosphere; i.e., it lifts the material out of the 
gravitational field into the dust-forming zone. With increasing 
strength of the pulsation amplitude Am, the additional input of 
energy and momentum leads to an enhanced density in the dust- 
forming region and supports further dust formation rather than 
accelerating the wind. Consequently, the stellar wind becomes 



far more massive, and the mass loss rate increases from (M) = 
l.OxlO^Moyr 1 for Am = lkms" 1 to(M) = 7.4xl0 -5 Moyr" 1 
for Am = 9kms~' at nearly constant velocity. As the wind be- 
comes more intense, the periodic depletion in dust-forming ma- 
terial also increases, and it takes more time to re-enrich the nu- 
cleation zone. The result is a longer shell period as seen in Table 

13 

Table 3. Parameter study for various excitation strength Am with 



fixed period P = 


376 d. 










Au 




/cDS 




<M> 


<Koo> 


<Pd/Pg) 


(kms- 1 ) 




(10- 3 d" 


') 


(M Q yr" 1 ) 


(kms- 1 ) 




1.0 


1.33 


2.66 


3.99 


1.0 x 10" 5 


26.4 


4.1 x 10- 3 


2.0 


0.53 


1.06 


1.33 


3.8 x 10~ 5 


29.9 


4.3 x 10 -3 


3.0 


0.46 


0.88 


1.33 


2.5 x 10 -5 


25.9 


3.3 x 10~ 3 


4.0 


0.47 


0.90 


1.35 


8.5 x 10- 6 


25.6 


3.2 x 10 -3 


5.0 


0.47 


0.93 


1.33 


2.1 x 10- 5 


28.2 


3.7 x 10~ 3 


6.0 


0.46 


0.79 


1.30 


4.1 x 10- 5 


30.0 


4.0 x 10~ 3 


7.0 


0.43 


0.93 


1.29 


1.9 x 10~ 5 


28.4 


3.7 x 10 -3 


8.0 


0.45 


0.86 


1.33 


1.2 x 10~ 5 


28.4 


3.5 x 10~ 3 


9.0 


0.44 


0.87 


1.69 


7.4 x 10~ 5 


30.6 


3.5 x 10" 3 



Values in boldface indicate the most dominant mode, whereas sin- 
gle and double underlined values refer to the excitation frequency 
and its first harmonic, respectively. 



Figure|7]depicts the spectrum and maps of the reference CDS 
excited with period P = 376 d and amplitude Am = 9kms~ 1 . 
With increasing amplitude the CDS behaviour changes from 
multiperiodic to irregular, as can be seen in the stroboscopic 
maps (cf. Fig. [6] lowermost panels, right and Fig. [7] right). The 
double periodicity breaks up. As can be seen in the power spec- 
trum, the timescale of the stellar pulsation (excitation modes) 
imprints their dynamics onto the envelope close to the star. With 
onset of the dust formation at r ~ 2.5 R* a continuous power 
spectrum appears rather than distinct modes in terms of the ra- 
dial excitation frequency. 

5. Summary and outlook 

In this article we have investigated the dynamics of envelopes 
around standard and low-luminous, carbon-rich AGB stars. We 
analysed the numerical solutions of the coupled equation sys- 
tem of hydrodynamics, thermodynamics, dust nucleation and 
growth, and radiative transfer by means of power spectra and 
stroboscopic maps of quantities, such as radial outflow velocity, 
temperature, and density. 

As shown for the case of CDSs around highly luminous 
AGBs, the presence of dust influences the dynamical and ener- 
getic behaviour of the entire system. Firstly, it amplifies the mo- 
mentum coupling between matter and radiation field, and sec- 
ondly, it enforces an additional dynamical behaviour (exterior 
A--mechanism) besides the pulsation to the dust envelope. The 
consequence of the coupling is the development of an envelope 
dynamics, which is in fact excited by the momentum and energy 
input by the stellar pulsation, but is operating on the timescales 
of the physical processes that determining the envelope. These 
dynamics can be characterised by a set of autonomous frequen- 
cies that are not identical to the frequency of the pulsation. 
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Fig. 7. Power spectrum and stroboscopic maps for a CDS excited with a period P = 376 d and strength Au — 9 km s 1 . The plot 
parameters for the maps are the same as in Fig. [6] 



For CDSs around low and standard luminous AGBs, such 
distinct frequencies are not obvious. To study the shell eigen- 
dynamics, the system has to be supplied with energy and mo- 
mentum. This we have done by introducing an additional force, 
regarded stochastic, e.g. provided by convection in the giant's 
atmosphere. The goal is to supply enough energy to introduce a 
CDS generating mass loss, but also to minimise the influence on 
the CDS's eigendynamics. 

When applying this method with different stochastic noise 
intensities to a reference CDS, the power spectra of the radial 
velocity show up as a set of eigenmodes of about 10 3 d, which 
represents the dynamical timescale of dust formation and el- 
ement enhancement in the dust formation zone. To study the 
interaction between the envelopes's eigendynamics and stellar 
pulsation, a parameter study of various stellar excitation peri- 
ods and strengths was carried out. Depending on the stimula- 
tion, the system reacts periodically, multiperiodically, or irregu- 
larly. For short excitation periods, the considered power spectra 
suggest that the CDS models tend to be semi-regular or even 
chaotic. The timescale of the excitation is much smaller than 
the timescale of the various physical processes (in particular the 
timescale of dust formation) in the envelope, so that the shell 
is not able to follow the excitation. For long periods the system 
is dominated by the external excitation. In contrast to envelopes 
around high luminous AGBs, no eigenmode-dominated domain 
could be found. In a closer examination of excitation periods 
from the irregular domain, we found the same resonance phe- 
nomenon (iP - jPcDs) as in the case of envelopes around high- 
luminous LPVs and Miras. However, it turns out that the integer 
j always equal 1 for standard luminous AGBs. 

Combined with the analysis method of shells around highly 
luminous AGBs already presented in a previous article, we are 
now able to start a systematic stellar parameter study. In order 
to compare the results with observations, this planned study will 
also analyse synthetic lightcurves and spectra. 
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